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ABSTRACT 


The work done in this investigation is part of a continuing program 
designed to develop probabilistic strength and life models for filament 
composite materials. In this particular experiment, groups of single fibers 
from two chemically identical different production spools of graphite were 
loaded to identical sustained tensile loads in order to determine their life 
Statistics while under stress. While previous work has produced models 
of composite reliability as a function of strength, this work develops a 
model which will be necessary to predict the reliability of composite fibers 
as a function of time. Parametric and non-parametric methods were used 
to quantify the relationship between fiber reliability as a function of stress 


and as a function of time. 
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I. INTRODUCTION 


While still in its infant stages, the use of composite materials in 
manufacturing is steadily increasing, particularly in high reliability, low 
accessibility environments such as aircraft fuselages and satellite fuel tanks. 
High strength to weight ratios, combined with minimal maintenance 
requirements and extreme resistance to corrosion make composites ideally 
suited for structures previously built exclusively of steel or aluminum. 
However, because composite usage has such a short history, there is no 
significant data base which would allow accurate reliability predictions. 

Through an ongoing program of research and experimentation, The 
Advanced Composites Laboratory at the Naval Postgraduate School at 
Monterey, CA, is defining a probabilistic model which 1s based on the 
failure mechanism of composites in eusian. Small scale designs which 
exploit fiber strength under tension, have accurately verified the Local 
Load Sharing model. More recent work is concentrating on the application 
of this model to large scale structures. With this in mind, it is important to 
understand that a structure must not only be predictable (reliable) in 
strength, but also in life. Because of their non-homogeneity, the study of 
composite materials must be further broken down to the study of not only 
individual constituents, but also of their inter-relationships with each other. 
For our purposes here, all work 1s done exclusively with graphite/ epoxy 
composite.; therefore, a proper understanding of the fiber, the matrix, and 
the interaction between them, is essential to formulating an accurate model 


of composite reliability. 


The primary objective of this project was to examine a large sample 
space of graphite fibers under different load conditions, and formulate a 
model of their reliability in life. The fiber life results would then be 
applied to composite life experiments now on-going designed to predict the 
life of composite structures. With data from these two projects in hand, 
one could then ostensibly isolate the solution to the fiber/ matrix interface 


problem and produce a realistic life prediction model. 
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Il. BACKGROUND 


A. STRENGTH DETERMINATION OF FIBER REINFORCED 
COMPOSITES 


Fiber Reinforced Composite materials generally consist of high 
strength fibers bonded together in a ductile matrix. In these materials the 
fibers carry essentially all of the stresses imposed on the structure while 
the matrix serves as the primary vehicle for transferring load between 
adjacent fibers, known as the Local Load Sharing model. In this model, 
the local load carried by a bundle of fibers is carried equally until the 
weak fibers begin to rupture. When a fiber breaks, its load 1s transferred 
to the adjacent fibers through the matrix. Rosen has modeled the load 
transfer by the existence of an ineffective length at the end of each 
(failed) fiber [Ref. 1]. In this region, the longitudinal stress (6) diminishes 
to zero while shear stresses (T) rise to a maximum. As these ineffective 
lengths accumulate, loads continuously transfer from the weak (failed) 
fibers to the strong ones. If enough fiber failure sites are in close 
proximity to each other, the matrix begins to crack or de-bond from the 
fibers. When this occurs, the matrix can no longer transfer the load and 
the sample fails catastrophically. Thus composite strength can be 
modeled as a function of fiber strength, the ability of the matrix to transfer 
load, and the interface between the fiber and matnix. The mathematical 
model which predicts the probability of failure of the composite (in 
tension) from the statistics of fiber strength and the effective length has 


been established by Harlow and Phoenix [Ref 2 and 3]. 


These three quanta are recognized as building blocks in composite 
strength determination. Matrix optimization is guided by this model. From 
a given assortment of matrix materials, one may determine which has the 
most suitable characteristics for stiffness and transferring shear stresses. 
Fiber strength statistics are determined through loading to destruction 
measurement. However, it is not feasible to directly measure and quantify 
the molecular level interface region. This can be accomplished by 
measuring the strength of a composite of known fiber and matrix. the 
fiber strength statistics and composite strength statistics can then be used 


simultaneously to solve for the ineffective length parameter. 


B. LIFE DETERMINATION OF FIBER REINFORCED 
COMPOSITES 


Decades of experience with homogeneous materials support the 
correlation between strength and life: that 1s, the larger the factor of 
safety for a given structure, the longer it can be expected to last. 
However, Coleman's theory of time dependence in the mechanical 
breakdown of fibers [Ref. 4] implies that the strength/life relation in load 
bearing fibers may not be so straightforward. Using the “weakest link 
hypothesis" first developed by Pierce in 1926, Coleman concluded that 
fiber life was instead a function of the statistical distribution of these 
“weak links" in tme dependent strength. The time dependent strength of 
the fiber, termed breaking kinetics, 1s the extreme value of fiber life under 


a given load history. This is fundamental to the modeling of composite life. 


Only when fiber life statistics are available can the model be applied to 
composite structures. Hence, the purpose for this work. 

The theory proposed is, the strength of a specific fiber can be related 
to the life of that same fiber. The dichotomy of course is, once a fiber has 
been broken in the strength test, another molecularly identical sample is 
not available for the life test. The alternative is to obtain a set of 
statistically identical fibers and test a portion of them in strength and 
another portion in life. 

This is being accomplished by using two different spools of AS-4 
graphite fiber, AS4-008 and AS4-019. A statistically meaningful number 
of samples have been taken from these two spools to estimate the 
respective parameters of the Weibull Distribution. The estimated 
parameters of these two different spools are different, suggesting the 


intrinsic strength of these two spools are from two different populations. 


Il. FILAMENT LIFE TESTING METHOD 


A. TEST SPECIMENS 

Test samples were single strands of Hercules Magnamite AS4 cut to 
two inch gauge lengths. Each sample consisted of one fiber filament 
mounted on a paper cardboard with a cutout for handling. The cardboard 
mount was severed with a hot wire prior to load application. Specific 
handling and mounting procedures are described in Appendix A and in 
Reference 5. Once mounted, the fiber diameters were measured using a 
laser diffraction device also following the procedures described in 


Reference 5. 


B. TEST EQUIPMENT DESCRIPTION 

The single filament samples were arranged in groups of 64 at four 
different sustained load levels by dead weights. The elapsed time 
between the startup of load and ending of stress-rupture at each station 
was triggered by an infra-red light switch. Test equipment is described in 
Appendix B. A polymethyl methacrylate (PMMA) cover was used to 
protect the samples from ambient air disturbances and ultraviolet (UV) 
radiation. It is known that UV light degrades the strength of some fibers; 
therefore, PMMA was chosen as a shield and is analyzed in Appendix C. 
Once loaded, the fibers were continuously monitored by a dedicated PC 
(personal) class computer. Software programs are reproduced in 


Appendix D. 


C. TESTING PROCEDURE 
1. Strength Test 

Fiber strength data were acquired previously and are covered in 

References 5 through 7. 
2. Life Test 

During the diameter measurement process, every fifth fiber was 
tested to failure in strength using a calibrated load cell. This provided the 
link between the two groups needed to verify that the samples tested in 
life were statistically identical to those tested previously in strength. The 
actual loading of the samples on the board was carefully controlled and, 
as much as possible, were all done by the same operator. The loading 
procedure is also described in Appendix A. 

There were four levels of load set up on the board, 14.17 g. 
eeve! 1), 3.15 s. (Level 2), and 12.20 g. (Level 3) and 10.96 g. (Level 
4). Each load level had 124 stations, 64 for each spool. This allowed a 
total of 512 life samples to be measured at sustained loads. Level 4 was 
planned but postponed. Instead those stations were also loaded to 14.17 
g. to gain additional measurements at that level. The sample board 
alternated rows with 008 and 019 fiber to give 4 rows of each. After 
each fiber was loaded, the date and time of the load were recorded in the 
lab notebook and also in a software spreadsheet. If a fiber failed on 
loading, its failure stress was recorded. After loading, each station was 
actively monitored by an AT computer which recorded stress failures as 
they occurred. Exact failure information was periodically retrieved and 


recorded in the notebook and spreadsheet. The spreadsheet 


automatically computed realized data for stress ruptures, early censor 
(Type Iv) stress ruptures, fail on load stresses, and existing lifetimes for 
each fiber type at each load level. It would then arrange the desired data 


in a form easily transported to other programs for analysis or presentation. 


IV. RESULTS 


A. TEST DATA 
1. Overview 

The data presented herein will be used to identify a suitable 
strength-life model for fiber filaments. A probabilistic formulation of the 
strength-life model is a joint distribution of the random variables, stress 
(applied), o,and elapsed time, t, denoted f(0,t). The parameters of this 
joint distribution are determined by three classes of realized data. The 
first is data realized in strength over negligible time. The second is data 
realized in life for a given stress history, in this experiment, constant, 
sustained stress. The third is data realized in strength given a survival in 
time history. The first two cases, that is, strength and life data, will be 
presented in this investigation. To examine variability in production, two 
spools of AS-4 graphite fibers were treated as two distinct populations, 
designated as spool! 008 and 019. The division of fibers from these two 
spools is traced out in Figures | and 2. The entire population of either 
spool was denoted as {D}. This population is divided into strength test 
samples yielding data, {D,}, and life test samples yielding data, {D; }. The 
respective sub-division of these two data sets 1s described in the following 


two sections. 


2. Strength Data 
The strength data set comes from work done previously using the 
same 008 and 019 fiber populations [Refs. 5 through 7]. The sources of 


this data are summarized in Table 1. 


TABLE 1. SOURCES OF FIBER STRENGTH DATA 


DATA NAME REF. NO. 
SE 

{D,;} Dr. Edward M. Wu 6 
Mr. Glen Nypiak 

Ores Lt. David Bell, USN 5 

(De LCDR Carl Englebert, 7 
USN 

{Dgq } Mr. Jim Nageotte 7 

3. Life Data 


Within the life data set, {D; }, not all of the samples were realized 


in life. Some of the samples were tested to assure the life samples were 
Statistically the same as the strength samples and are denoted as {D;,} 
(Figures | and 2). Some samples failed before they reached the target 
sustained load level. Thus, this was part of the life test subset but was 


realized in strength and is denoted as {D+ ;,}, {D+ 2, }, and {D,z,}. Based 


on kinetic theory [Ref. 2], it is generally believed that the load supported 


by a fiber increases as the loading rate increases; therefore, failures in 
strength would be hypothetically eliminated if the fibers could be loaded 
instantaneously. While this not physically possible, the data from samples 
which failed on loading can be accounted for by treating the life samples 
conditionally to strength greater than the sustained load. Data obtained 
from this experiment is summarized in Tables 2 through 4. Ruptures 
occurring on loading, prior to the desired load are identified as "Fail on 
Load". Stress ruptures occurring in time but caused by other than 
intrinsic stimulus internal to the sample are “Type Iv’ censored. This data 
is Statistically treated as greater than or equal to the realized life. The 
censorings are chance events as denoted by the variable time at censor. 
Stress ruptures in time are "Exact" and fibers still under load at the end of 
the experiment are "Type I" censored. This data is statistically treated as 


greater than or equal to the time the data was interpreted. 


fame 2 RESULTS OP EOAD LEVEL | 


Type of Failure AS4-008 AS4-019 
Fail on Load {D; j,} 53 on 
Type Iv Censor {D, 1,,} 9 7 
Exact {Dj 44 } 15 A) 


Type I Censor {D} 4, } 116 110 


TABLE 34 RESULTS OF LOAD EE Ere 


Type of Failure AS4-008 AS4-019 
Fail on Load { D9, } 26 13 
Type iv Censor 4 ] 
{Dray} 

Exact {D5} } 9 6 
Type I Censor {D; 9, } 62 63 


TABLE 4.. RESULTS OF LOAD BEY rire) 


Type of Failure AS4-008 AS4-019 
Fail on Load {D,3z, } 10 1] 
Type Iv Censor ] 3 
{Dizy} 

Exact {D;3; } 3 | 
Type ) Censor {Diz } 62 63 


B. DATA ANALYSIS 
1. Analysis as a Function of Stress 
a. Overview 
Before proceeding with an analysis of life based on strength, 


it must first be confirmed that the two sets of data, {D, } and {D, }, are 


indeed statistically from the same population. As can be seen in Figures | 


and 2, each spool has yielded two sets of samples. The first set, {D, }, 
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was tested entirely in strength. Within the second set, {D; }, every fifth 
fiber was tested in strength in order to link the two sample sets {D,} and 
{D; }. 
b. Non-Parametric Interpretation 
Non-parametric data analysis was conducted by plotting the 
cumulative failures in a weakest link (Weibull Probability) coordinate. 
This required the data to be ordered and a rank assigned. For this 


experiment, Expected Rank was used and is defined: the if realization of 


the ordered x; 1S 





Ol ase 
N+] 


and x; is the ordered data (from weakest to strongest) in a set of N 
samples. An ambiguity arises when one considers the inclusion of 
censored data. Because the point 1s censored, its exact strength, and 
therefore its exact rank, cannot be precisely determined. Furthermore, 
because the data is ranked, this uncertainty perturbs the rank of each 
subsequent data point greater than the censored data. For this reason, 
censored data was omitted from the graphical representation. However, 
censored data did not affect the Likelihood calculations and was therefore 
included in the parametric interpretation. 

Once the data has been ranked, it is linearized by using the 
weakest link transformation. This method linearizes the ranked data and 
allows a visual observation of the lower tail characteristics. The weakest 
link transformation 1s defined as: 


F*(o) = In(-In(1-F(o))) 


where F(O) is the empirical rank of the data in question. Recalling that 
each fiber is modeled as a chain of independent links, the strength of a 
given length of fiber must be the same as the strength of the weakest link 
in that chain, hence F*(o) can be further defined as the probability of 
failure of the weakest link of a given gauge length of fiber. 

Strength data, {D,}, from these two spools (008 and 019) is 
plotted in Figures 3 and 4. This set of fibers had previously been 


confirmed to have come from the same set of like materials and gauge 


lengths (Ref. 9). Figures 5 and 6 depict the F* plot of the set Dy, overlaid 
with the set Ds. The close proximity of the two sets of data indicates they 
too are of the same fiber and gauge length. A cursory run-test on the 
ranked set of merged data shown in Figures 7 and 8 shows acceptable 
levels of statistical clumping and completes the non-parametric 
identicalness test of the two sets. 

c. Parametric Interpretation 

A parametric analysis of the same data was also conducted 


by using the Weibull distribution for both strength and life: 


a 


x 
F(x) =1-exp;-+ — 
B 


where a and B are the shape and location parameters respectively. 
Given a data set, these parameters may be obtained by using the 
Maximum Likelihood Estimator (MLE) derived in Appendices E and F. 
The attribute of this method is there is no requirement to order the data 


which eliminates the process of ranking which allows the inclusion of any 


censored data points. Also associated with the MLE is the concept of 
confidence interval. The confidence region of the parameters a and B 
can be obtained from the likelihood contour plots. The confidence region 
is the confidence, stated as a percentage, that a selected o and B are the 
actual @ and 8. Parameters chosen at a very low confidence interval 
have a low probability of matching the actual parameters, but when they 
do, they are highly accurate. A very high confidence interval, while 
allowing more room for error, is not very precise and is of little practical 
value. 

Figures 9 and 10 show a three dimensional representation of 
the Likelihood for the D, data set along with various confidence region 
contours. The 5 % and 1 % confidence region contours are shown in 
Figures 1] and 12. The shape parameter, a, and location parameter, B, 
for this data set can be determined from Figures 3 and 4 graphically by 
the slope of the line of data points and the value of the X axis at F* = 0 
respectively. Alternately, the estimate can also be based on the calculus 
MLE derived in Appendix E. 

In order to show that {D, } 1s identical to {D;,}, the likelihood 
ratio can be used. The likelihood ratio is graphically represented by the 
volume of the intersection of the two likelihoods at a given confidence 
interval, divided by the volume of the same confidence interval of the 


accepted baseline data set, {D, } or, symbolically 


V(D,)f )V(D..) 
V(D,) 


: V(D,) and V(D;,) are the volumes of the likelihood contours 
within the desired confidence interval. Figures 13 and 14 show overlay 
plots of the Likelihood contours of N, and N;, at 5% and | % confidence 
intervals. By examining these figures, it can be determined that the 
likelihood ratios for 008 and 019 at the 5 % confidence interval are 
estimated to be .75 and .70 respectively, which is deemed to be 
acceptable for this confidence interval. 

The two sets of data can now be merged and the tested in 
life fibers are confirmed to be of the same population as the tested in 
strength fibers. As a last point of confirmation, Figures 15 and 16 overlay 
this merged data set {D,+D;,} with the original set,{D,} showing they 
are essentially identical, as do the contour overlay plots shown in Figures 
17 and 18. 

A distinction can now be made concerning those fibers 


which failed on load ({ D; 4,}, {D+ 2, }. {D+3, }). Although intended to be 


measured in time, these fibers actually failed in strength prior to reaching 


the target loads, 6,, 67, and Gz respectively; therefore, this set of data can 


be considered as realized in strength to be included in the merged strength 
data set {D,+D;,}. The Expected Rank of the exact failed in strength 
data was presented as before with the exception that, the samples which 
reached the target load in life test were considered Type I strength 
censored at the stress rupture load, This is the equivalent to stating its 
strength would be greater than the stress rupture load. As with a// non- 
parametric analyses, these censored points were not presented in the F* 


domain because of the ambiguity in rank. By treating censored data in 
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this manner, the graphical representation of the middle and upper tails 
remains unaffected, while the resolution of the lower tail is improved. Fail 
on loading data, {D;4,}, {Dio,}, and {D,3,}, is presented in Figures 19 
through 24 and can be analyzed in the same manner as the set {D;,} with 
similar conclusions being drawn. Thus, the set of fibers tested in time can 
now be considered statistically identical to the set of fibers tested in 
stress. 
The final step in strength analysis was to analyze the entire 
population of strength data and determine a shape (a) and location (B) 
parameter for the Weibull Probability of Failure Function which would be 
necessary to determine the Likelihood in Life. Results of non-parametric 
analysis of the set {D} in strength are shown in Figures 25 and 26. 
Because of the data censored above 12.20, 13.15, and 14.17 grams, a 
value for & and B can only be estimated from the set {D, + D;,}; 
therefore, a parametric analysis in the form of the Likelihood contours of 
the set {D} 1s shown in Figures 27 and 28. 
2. Analysis as a Function of Time 
Since the life of a fiber is dependent on the fiber’s surviving the 
loading process, ruptures during loading must be considered in 
determining the Likelihood of the fiber in life. The reliability of the fiber in 
strength will be necessary to determine its reliability in life. Therefore, 
equation (1) in Appendix E must be modified to include the fibers 


reliability in strength: 


hy! 





normond-(55) | 


and the Likelihood in the time domain becomes: 





where 6, 1s the target stress and k is the number of left censored fibers. 
By determining estimates of && and Bo from Figure 27 for Spool 008 and 
evaluating the Likelihood contour over a set range of a and B,, the 
Likelihoods for each of the three load levels can be plotted and are shown 
in Figures 29 through 31 which were developed using an estimate of 5.0 
for Go and 18.0 for Bg. A similar procedure was followed for Spool 019. 
Although not investigated during this analysis, an alternative 
method of calculating the Likelihood would be to evaluate the fail on 
loading data as censored in time rather than stress. In order to do this, 


One must first assign a time, tg, at which the sample is considered to have 


attained the target stress. All samples which fail on loading can be 


considered to have a life less than or equal to tg, or: 
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Under this interpretation, the Likelihood can be calculated using: 


k 
if | m mo Ge Bet so 
F(t) = - ens - (2) | Cay xS' expt B, oF] 
t i=] 


m RQ. . xX . 
Fi owl -(E] jew-(32) 


3. Construction of the Strength-Life Model 


(n—m) 


For a life distribution given a sustained stress level, o;, the Weibull 


Fitlo=o)=1- =P ce 


Using this model, the Likelihood contours for the life data realized up to 


model used was: 


the current time were plotted and are shown in Figures 29 through 31.The 
flatness of the contours is the result of an extremely broad range of B 
which 1s caused by the elapsed time (at the time of data interpretation) 
being much smaller than the underlying location parameter of life, B. 
Therefore, only very tentative estimates of B are possible at this time. 


The shape and location parameters are visually estimated to be: 


We 


A 


Stress Level, o; 8k o 


14.17 g. 15 1.0E+11 
lige pal lis ex is) 1.0E+12 
122 0re. 1S 5.0E+13 


The estimated life at different levels of sustained stress is shown 
in Figure 32 evaluated at three levels of probability of failure, 0.001, 0.01, 
and 0.1. The figure can be considered as a graphical representation of the 
Strength-Life Model for AS4 graphite under sustained stress. The 
expected life of the fiber, evaluated at the desired probability of failure at a 


given stress level, can be obtained directly from the model. 


V. CONCLUSIONS AND RECOMMENDATIONS 


The Strength-Life Model, once verified, can be a useful and important 
tool for predicting fiber reliability using strength statistics. In design, it 
allows accurate determination of maximum design stress for an expected 
life. In maintenance and repair, it facilitates the calculation of the 
remaining life of an existing structure. It also predicts stress levels at 
which an existing structure may be operated that can increase or 
decrease its remaining life. 

The following items are recommended for follow-on research:. 

1. Apphcation and verification of the model using larger scale 
composite samples. 

2. Mathematical formulation to predict life based solely on strength. 

3. The life test be continued to improve the estimation of the 


parameters of the model. 
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APPENDIX A 


PROCEDURES. 


Fiber Handling. Fiber handling procedures are specific and are in 
effect continuously. 

1. Unless required for the immediate task at hand, all fibers and 
filaments shall be stored in dark, secure containers. 

2. Fibers and filaments are fragile. Any fiber or filament that is 
dropped or otherwise subjected to extreme stress, shall be discarded. 

3. Whenever possible, all samples should be handled by one end 
only so as not to exert any tension on the sample. 


4. Use of fluorescent lights shall be kept to a minimum. 


Filament Mounting. Filaments were prepared according to the 


following guidelines: 


1. Strip off a random length of fiber from the spool, sever, and 
discard. 

2. Cut the next 12 - 15 cm. of fiber from the spool and tape one end 
to the glass plate. 

3. Pour a small puddle of ethyl alcohol on the glass plate to float the 
free end of the bundle. The filaments will begin to float apart. Gently 


fan out the filaments floating in the alcohol and let stand until all the 


i) 
NO 


alcohol has evaporated. 

4. Write the correct serial numbers on the carriers to be loaded. 

5. Using a small piece of transparent tape on the end of a modelling 
knife, lift the free end of one filament and slide an empty carrier under 
fe 

6. While maintaining filament alignment with the two small holes in 
the carrier, lower the raised filament end and tape it to the glass 
plate. 

7. Use two small pieces of tape to secure the filament to the 
Carrier. 

8. After the filament is secured between the two pieces of tape, 
carefully cut the filament to free it and the carrier from the glass plate 
and place the filament/carrier assembly on the cold curing plate. 

9. Repeat steps 5-8 eight more times. After the ninth filament has 


been mounted, discard the remains of the bundle. 


NOTE 
In order to maintain standard curing times and temperatures, do 
not activate the heat element in the hot plate until the plate 1s full and 
do not remove the samples until the plate has been turned off and has 


cooled. 
10. Activate the hot plate with a setting of 200° F. 


11. Thoroughly mix a small amount of epoxy (Devcon 2 Ton). 


12. Reduce hot plate temperature to 150° F. 
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13. Place one dot of epoxy at each end of the elongated hole in 
each carrier to permanently secure both ends of the filaments to their 
respective carriers. Also, place a dot on the outboard ends of the two 
small holes to re-enforce them. 

14. When the last filament has been epoxied, turn off the hot plate 
and let cool. 


15. Store mounted filaments in designated containers. 


Filament Loading. 


Loader Start up 


1. Turn on power supplies for load cell and elevator. 


ID 


Turn on HP-85 and attached Data Acquisition Unit. 

3. Insert program tape into HP-85 and type LOAD "QUIK4". 

4. Type RUN. 

5. Select "k2" (INPUT) key and enter the date. 

6. With load cell assembled and unladen, select "k!1" (ADJ B) key. 
This sets the bias to 0. To check, select "k3" (WEIGH). If a weight 
within .003 g of 0.0g is not obtained, re-select ADJ B until satisfactory. 


Accurate readings should also be verified using calibrated weights. 


Filament Loading 


I. Remove Plexiglass cover from sample rack. 


24 


2. Without disturbing adjacent filaments, lower stabilizer bolts as far as 
possible on stations to be loaded. 
3. Adjust the spring so that the flag will operate the optical trigger 
without rubbing against it. 
4. Remove the weight from the station to be loaded. 
5. Using tweezers, grasp the desired filament by one end of it's carrier 
and hang it on the hook attached to the flag. 
6. Using tweezers, hang the weight on the bottom hole of the carrier. 
7. Verify the flag does not rub against the trigger. 
8. Lower the load cell elevator far enough that it will fit onto the rack 
without disturbing the hung filament. 
9. Without anything touching the top of the load cell, press “ADJ B" on 
the HP-85. 
10. When the screen has stopped scrolling, raise the elevator until the 
compliance spring is completely retracted and the paper carrier 1s 
carrying exactly zero load. 
11. Verify that the PC records the station as unloaded as the flag 
activates the optical trigger. 
12. Select "LOAD" on the HP-85 and answer the cues: 

“ENTER THE SAMPLE #° 

"ENTER THE STATION #° 
fe, Select "CONT" 
14. Verify the weight of the vial shown by the computer is correct for 
the station being loaded. 


15. Ensure the direction switch on the elevator is in the "down" 
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position and the power switch “off”. 

16. Check the setting on the hot wire power supply and turn it on. 
The wire should glow a dull red, not bright. 

17. Without touching the filament or disturbing adjacent filaments, 
carefully burn through both sides of the paper carrier at the bottom of 
the slot. 

18. Turn off the hot wire and return it to it’s holder. 

19. Select "CONT" on the HP-85 and move the elevator power 
switch to "on". Monitor the elevator as it lowers to ensure it does not 
bottom out prior to being unloaded. 

20. Verify the PC records the station being loaded as the flag blocks 


the optical trigger. 
NOTE. If the filament breaks before it is completely loaded, stop 
the elevator, record the mass displayed by the HP-85, and start again 


from step 24. 


21. When the weight hangs free, the HP-85 will beep twice. When 


this happens, stop the elevator. After two minutes, the HP-85 will beep 


twice again, allowing you to continue. 

22. Lower the elevator until it can be withdrawn without disturbing the 
filament. 

23. Once the elevator has been withdrawn, carefully, without 
disturbing the filament, adjust the stabilizer screw so that it is 


centered on the vial without touching it and has a vial/stabilizer gap of 
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not more than | mm. 

24. Load a sheet of paper into the plotter and select "shift k2" (PLOT) 
on the HP-85, answering the cues as required. 

25. Record filament load information in Fiber Life Data Book and file 
load plot in the plot binder. 

26. Repeat from step 2 as necessary. 

27. When all desired stations have been loaded, verify all stabilizer 
bolts are positioned correctly and replace Plexiglass cover. 


28. Verify PC is reporting the correct status for the stations loaded. 


Us || 


APPENDIX B 


EQUIPMENT DESCRIPTION 


Data for the experiment was collected by suspending fiber samples on 
the individual stations (512 total) of a large sample board. Load and fail 
status was actively monitored by individual infrared switchesthrough 
ribbon wire which fed the data to an interface box communicating with a 


PC-AT. (Figure B-1). 


The sample board was made up of 512 stations arranged in 8 rows 
and 64 columns. Each station consisted of a mounting point, a high 
compliance spring (k=0.635 mm/g), an optical switch and trigger, a 
filament sample, a weight. and an adjustable stabilizer screw. The board 
was constructed of aluminum channel stock and fiberboard. When not 
being loaded, the stations were covered with a protective sheet of 1/4" 


PMMA. 


The trigger was a wire paper clip with the inner loop removed. The 
top half of this outer loop was then wrapped in black viny! electricians 
tape. The fiber sample, mounted in its paper carrier, was hung on this 
trigger which in turn hung from the,bottom of the spring. With the weight 
attached to the bottom of the fiber carrier, the black tape blocked the 
optical switch. When the filament broke, the weight dropped away while 


the spring retracted the trigger allowing light to pass through the paper 
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clip thus closing the optical switch, which was detected and recorded by 
ihe PC. 


The weights were plastic vials filled with lead shot and were loaded to 
3 weight levels: 12.20, 13.15, and 14.17 g. Two rows of AS-4 graphite, 
batch no. 008, and two rows of AS-4 batch no. 019 were loaded at 14.17 
g., one row of O08 and one row of 019 were loaded at 13.15, and one row 
of 008 and one row of 019 were loaded at 12.20 g. It was decided to 
double the number of high weighted samples in order to obtain more data 


points prior to censoring the experiment. 


Filaments were permanently mounted in paper carriers to facilitate a 
reliable and convenient method of storage and handling. Once assembled 
in itS proper station, the paper carrier was severed using a hot wire 
instead of scissors to reduce the chance of sending a shock through the 


filament. [Ref. 3] 


Each of the 512 optical switches was connected to a custom designed 
interface box. This box contained a power supply for the optical switches 
as well as a multiplexing circuit which fed status information to a 
MetraByte I/O card in the PC. The board was divided into 64 8&-bit 
words and software was written to cycle through each of the words and 
to report any changes in the value of the word. For example, if stations 1- 
8 are all loaded the PC reads a value of 0000 0000. If station 2 fails, the 
next time this byte is read, it would show a value of 0000 0010. The PC 


would then record the failure and continue monitoring the board. A back- 


I) 


up power supply was connected to the interface box and PC in case of a 
filament failure during a power outage. PC software is reproduced in 


Appendix D. 


Start times for each filament were also recorded by the PC when a bit 
change from | to 0 was detected. However, consistently accurate start 
times were highly dependent on a function of individual spring compliance, 
precise trigger construction , loader motor speed, and masses of the 
weights. In order to eliminate this variability, the loader software used to 
program the HP-85 also incorporated a ttming mechanism. As the load 
cell was lowered, the HP-85 monitored the weight on the cell. When zero 
weight was detected, a beep sounded and a two minute timer was 
started. The load cell was left under the sample until expiration of the two 
minute time period. Should a sample have failed during the first two 
minutes of life, its lifetime would be recorded to the nearest 1/10 second. 
After two minutes, such precise record keeping was not necessary and 


timing was turned over to the PC. 


The HP-85 was programmed to read the load cell approximately 5 
times per second during the loading sequence and could therefore provide 
accurate data on the loads as they were applied. After loading, this 
information was sent to a plotter and a record of each load made for 
analysis. If a filament ruptured any time during the load sequence, up to 
the expiration of the two minute time period, the HP-85 would record the 
load on the filament at the tme of failure or, if already loaded, the life of 


the filament up to the first two minutes. 
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Filaments were loaded by using a 150 g. capacity Sensotec load cell 
mounted to an electrically controlled, hydraulically actuated elevator. The 
electric motor, Western Gear model P5B24R3 incorporated a built in 
reduction gear and was run at 7.5 volts giving a start to finish load run of 
about 30 seconds yielding 170-220 data points as the load on the filament 
increased. The motor was mechanically connected to a disposable 5 cc 
hypodermic syringe which acted as the master cylinder. The slave 
cylinder was an identical syringe with the load cell mounted on a platform 
secured to the plunger. This arrangement was in turn mounted to a 
Plexiglass platform which provided a stable base during loading. The 
slave and master were connected with 3 feet of 1/8" ID vinyl tubing; the 


hydraulic fluid was tap water. 


3] 


AS4-019: 14.17 g. 


AS4-008: 
AS4-019: 


AS4-008: 
AS4-019: 


AS4-008: 
AS4-019: 





AS4-008: 


Interface 
Box 





Figure B-1. Diagram of Fiber Life Testing Apparatus. 


APPENDIX C 


ULTRAVIOLET (UV) LIGHT SHIELDING 


The life of certain fibers, (Kevlar for example) are known to be 
affected by exposure to UV radiation. To assure these life test 
experiments were performed in an inert environment, an investigation into 
the presence of UV light was conducted. 

The sample rack was assembled in a secure laboratory environment 
on a basement level with windows facing Northeast. The rack itself was 
mounted on the Southeast wall and therefore, could not receive any direct 
sunlight. Windows were covered by opaque commercial grade drapery. 
Room lighting was fluorescent overhead with plastic covers. The sample 
racks were covered with 1/4 inch Polymethyl methacrylate (PMMA) 
Plexiglas. Light from the two sources was filtered twice before reaching 
the samples. The indirect sunlight passed through standard grade soda 
lime window glass and then through the PMMA rack covers. Overhead 
lighting was filtered through a plastic optical diffuser and then through the 
PMMA rack cover. 

Tests for the intensity of 300-400 nm UV light were conducted using a 
100 Watt mercury lamp at a distance of 30 ft. from the spectrometer. All 
measurements were made at the surface of the sample material. Results 
of the spectrum analysis are shown in Figures C-1 through C-3. To 


explore materials for improved UV shielding, analysis was also done on 
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Polycarbonate (Lexan) and a UV filtering grade of PMMA (Acrylite OP- 
2). Results from these two materials are shown in Figures C-4 and C-S. 
The improvement in shielding appeared to be a small range of an 
additional 25 nm to that already provided by the PMMA and was 
therefore not deemed to be significant. 

Based on these measurements, the amount of UV light reaching the 


actual fiber samples was judged to be insufficient to affect life testing. 
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Figure C-1. Results of UV Transmittance test for Soda Lime Glass. 
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Figure C-2. Results of UV Transmittance test for PMMA. 
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Figure C-3. Results of UV Transmittance test for Installed Flourescent 


Lighting Cover. 
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Figure C-4. Results of UV Transmittance test for Polycarbonate. 


38 


(WU) HLONZTSAYM 


00S OSV OO0V OGE OOE OSe 


ee ee ee Ae 3 eee oP 0° 


ees i ye ft hl 


(Z-dO aiAoy) VWI Buvari4 AQ :eoueyiusuesy 


ep) 


T 


Te) 


(0 


_ 


08) 


Oo) 


S$1UN09 


Figure C-5. Results of UV Transmittance test for UV Filtering PMMA. 
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APPENDIX D 


MAIN PROGRAM: 


'k* Filamod7 ** 
' Jim Nageotte Nov 13, 1990 


Start: 
DIM newrd%(64), O1dRd%(64), row%(8), 
Status$(2), StarTime$(512) 


+ 


Dimension Variables 


+ 


row%(1) = 1: row%(2) = 2: row%(3) = 4: * Declare constants 

row%(4) = 8: row%(5) = 16: row%(6) = 32: 

row%(7) = 64: row%(8) = 128 

a%¥ = 0 

OUT &H303, &H10 ' * Reset I/O Card 

TYPE Status2 * Declare file type 
OldRd AS INTEGER * for status file 

END TYPE 

DIM Stat AS Status2 

TYPE Times * Declare file type 
StarTime AS STRING * 16 * for Start times file 

END TYPE 


DIM Tim AS Times 


CLS 
PRINT “Ver <i1.13.90" 
GOSUB Gettime 


PRINT “Power up at "; TIMES; " on "; * Print start up message 
DATES * & current date on screen. 
PrintStat = INP(&H3BD) ' * Read printer status. 

* 223 indicates that the 

IF PrintStat = 223 THEN ' * printer is ready. 
LPRINT “Power up at "; TIMES; 
won": DATES * Print start up message 
LPRINT : LPRINT * only if printer is ready. 

END IF 

OPEN "Status.Dat" FOR RANDOM * Open Status File for 

AS #1 LEN = 2 : * random access. 

FOR idx = 1 TO 64 ! * Read previous status 
GET #1, idx, Stat * from disk file. 
OldRd%(idx) = Stat.OldRd 

NEXT idx 

CLOSE #1 * Close Status File. 

OPEN "Times.Dat" FOR RANDOM AS #2 * Open Start Times File 

LEN = 16 * and read start times 

FOR Pee) TO 512 * for all 512 stations. 

GELtz27° 13) Fim 
StarTime$(i) = Tim.StarTime 
NEXT i 
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CLOSE #2 


DO 


GOSUB Checksw 
FOR i = 1 TO 64 


IF newrd%(i) = O1dRd%&(i) THEN GOTO Same 
a0% = O1dRd%&(i): al% = newrd% (i) 
FOR j = 1 TO 8: jdx% =9 - j 

bO& = FIX(a0% / row% (jdx%) ) 

b1l%é = FIX(al% / row%(jdx) ) 

IF bO% = bl%¥ THEN GOTO Nextbit 

Station = (i - 1) * 8 + jdx% 

IF b1l% < 0 OR b1l%¥ > 1 THEN 
PRINT "Error in Status bl1%" 


BEEP 
END IF 
IF bl¥ = 1 THEN 

PRINT "Station #"; Station; 

"Loaded "; RdTime$; " on "; 

RdDates 

GOSUB Store 

PrintStat = INP(&H3BD) 

IF PrintStat = 223 THEN 
LPRINT "Station #"; 
Station; " Loaded "; 
LPRINT RdTimeS; " on "; 
RdDates$ 

END IF 

BEEP 

END IF 
IF b1% = O THEN 


GOSUB Elapsed 
PRINT: PRINT "Station #"; Station; 
"6"; StarTimesS (Station) ;"to";PRINT 


LEFT$ (RdTime$, 5); " ";RdDate$: PRINT 
" "; YEARS; "Years, "; DayT; "Days "; 
PRINT Min; " Minutes"; PRINT 
PrintStat = INP(&H3BD) 
IF PrintStat = 223 THEN 
LPRINT " ": LPRINT "Station #"; 
Station;" Ae 
StarTime$ (Station) ;" to "; 
LEFTS$(RdTime$, 5);" "; RdDate$ 
LPRINT " ";YEARS; "Years, "; DayT; 
LPRINT Min; "Minutes". 
LPRINT " " 
END IF 
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* + 


re 


® 
® 


++ + + + 


* + + + 


* Close Start times file. 


Read 8 bit status words. 
until something changes. 


Compares Status read from 
disk to that read from 
the switches. If they are 
not equal it looks more 
closely at the switches. 


0 or 1 are the only valid 
values for bl%. 


indicates "Loaded" 
Print loaded message to 
screen. 


Store new Start time to 
disk file. 

Check printer Status. 
If printer is ready 
then print loaded 
message to printer. 


"O" indicates failed. 
Find elapsed time. 
Print failure information 
to screen. 


Check printer Status. 
If printer is ready, 
then print failure 
information to printer. 


"Days " H 


END IF 


Nextbit: 


Same; 


a0% = a0% - bO& * rows(jdx%): al% = 


NEXT’ 3 


OldRd&(i) = newrd%(i) 

OPEN "Status.dat" FOR RANDOM AS 
#1 LEN = 2 

Stat.OldRd = O1dRd$%(i) 

PUT #1, i, Stat 

CLOSE 


NEXT i 


LOOP 
END 


SUBROUTINES: 


al% - bl& * row%(jdx$) 


* Update disk Status file. 


x End of loop that compares 
* old status to new and end 


* 


of main program. 


Check Switch compares the current status to the status read off the disk 
It then sets a flag and 
returns with the current time and date to start checking each of the 512 


file. 


stations and determine what type of change has taken place. 
the status in 8 bit words and in such a small loop, 


It loops until there is a change. 


By checking 
we keep the time 


required to check all switches and the response time to a minimun. 


Checksw: 


DO 
GOSUB Gettime 
at = 0: idx = 0: RdTimeS = TIMES 


RdDate$ = DATES 
FOR i = 1 TO 8: OUT &H301, row%(i) 
FOR 5 = 0 TO 7: OUT &H302, 5 
idx = idx +1 


newrdt(idx) = INP(&H300) 
IF newrd%(idx) <> OldRd% (idx) 
THEN at = 1 
NEXT 3 

NEXT i 

LOOP UNTIL at = 1 

BEEP 

RETURN 


+e + 


* + 


Get current time. 
Get current date. 
Select Row to read. 
Select word to read. 


Set Flag if Status has 
changed. 


Loop until Status has 
changed. 


Elapsed Time calculates the time between the Start time read from disk 
and the current time. It calls two other sub routines to-do this, Julian 
and Minutes. Julian is called twice, the first time it converts the 
Start time to a Julian date and the second time it converts the current 
date to a Julian date. The conversion is necessary to allow subtraction 
to be performed on the dates. Having arrived at the number of Years and 
Days before the fiber ruptured, we call the subroutine Minutes to 
subtract the times involved and return with the correct number of 
minutes. An adjustment to the number of days may be made based on the 
number of minutes. 


Elapsed: 
ConvDate$=RIGHTS (StarTimeS$ (Station) ,10) * Select Start Date. 
GOSUB Julian * Convert to julian date. 


StartYr = Year 

StartDy = Days 

ConvDateS$ = RdDate$ * Select Current Date. 
GOSUB Julian * Convert to julian date. 


StopYr = Year 
StopDy = Days 
YEARS = StopYr - StartyYr * Subtract years. 
DayT = StopDy - StartDy * Subtract days. 
IF DayT < 0 THEN 

DayT = DayT + 365 

YEARS = YEARS - 1 
END IF 
GOSUB Minutes * Calculate minutes. 
IF Min < 0 THEN 

Min = Min + 1440 

DayT = DayT - 1 


END IF 
RETURN 
Julian: 
Mon = VAL(LEFTS$(ConvDate$, 2) ) * Get month, day and year 
Day = VAL(MID$(ConvDate$, 4, 2)) * from the date string. 


Year = VAL(RIGHT$(ConvDate$, 2)) 
SELECT CASE Mon 


CASE 1 
Days = Day 
CASE 2 
Days = Day + 31 
CASE 3 
Days = Day + 59 
CASE 4 
Days = Day + 90 * Calculate total days. 
CASE 5 
Days = Day + 120 
CASE 6 
Days = Day + 151 
CASE 7 


Days = Day + 181 
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CASE 8 
Days = Day + 212 


CASE 9 
Days = Day + 243 
CASE 10 
Days = Day + 273 
CASE 11 
Days = Day + 304 
CASE ELSE 
Days = Day + 334 
END SELECT 
RETURN 


Store simply opens the Start Time file on disk, a random access file, 
and stores only the new time without reading or re-writing the 
unaffected times. The file is then closed and the routine exited. 


Store: 
OPEN "Times.Dat" FOR RANDOM AS #2 LEN = 16 
Tim.StarTime = LEFT$(RdTime$, 5) + " " + RdDate$ 
PUT #2, Station, Tim 
CLOSE #2 
StarTime$(Station) = Tim.StarTime 
RETURN 


Minutes recovers the values of Minutes and Hours from the time strings 
and converts them to Minutes only for ease of mathematical operations 
in the time keeping. | 


Minutes: 
HrSt$ = LEFTS$(StarTime$(Station), 2) 
Min$ = MID$(StarTime$(Station), 4, 2) 
Hrs = VAL(HrSt$) 
MinSt = VAL(Min$) 
MinS = Hrs * 60 + MinSt 
Hrf$ = LEFT$(RdTimeS, 2) 
MinF$ = MIDS$(RdTimeS$, 4, 2) 
Hrf = VAL(Hrf$) 
MinFn = VAL(MinF$) 
MinF = Hrf * 60 + MinFn 
Min = MinF - Mins 
RETURN 


Gettime reads I/O ports to read the time directly from the system clock. 
It then converts these values to a string that resembles the "TIMES" 
function. The program then uses this string called "Rtime$" for Real 
Time String. 
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Gettime: 
Sec$ = LTRIM(STR$(INP(&H2CO) AND &HF) ) 
Secs$ = LTRIM(STRS$(INP(&H2Cl) AND &HF)) 
Minute$ = LTRIM(STR$(INP(&H2C2) AND &HF) ) 
Minutes$ = LTRIM(STRS$(INP(&H2C3) AND &HF)) 
Hour$ = LTRIM(STRS$(INP(&H2C4) AND &HF)) 
Hours$ = LTRIM(STRS$(INP(&H2C5) AND &HF)) 

Rtimes = "00:00:00" 


MID$S(Rtime$, 1) = Hours$ 
MIDS(Rtime$, 2) = Hours 
MIDS$(Rtime$, 4) = Minutes$ 
MIDS (Rtime$, 5) = Minute$ 
MID$(Rtime$, 7) = Secs$ 
MID$(Rtime$, 8) = Sec$ 


RETURN 
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Read Data Port 






Determine which 
station has changed 


If loaded, display 
Sta. No. and start tume/ date. 
Record start time/ date 
on floppy drive. 


If not loaded, display Sta. No. 
and start/ fail time and date. 
Record Sta. No. and fail 
time/ date on floppy drive. 





Figure D-1. Flowchart for Life Test Monitoring Program. 
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APPENDIX E 


I LIKELIHOOD DERIVATION FOR EXACT DATA 
For Exact Data: 


eee = xe, B) = d/ax F(x., o, 8) 


-F(x,0,8) =1 -ex p- (+) | 


Taking the derivative with respect to x: 
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The Likelihood for all data, x;, = 


n 
i)" 


me 84 
= op 
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() 


(2) 


fet A = In = In [(x,,c,6) 


- Inf a" ae. exp4 - pax: b] 


2 (a-D -a a 
=nIna-nalngp+ Sin x. exp - B ze} 


=n Ina —nalnB+(a- Din x f a inexp rs perth 


=a). roe (3) 


n 
=nIna-nalnB+(a—- Din Nn 2x. 
i= 


To find the maximum Likelihood, (L,,,, ), set the derivative of A with 


respect to « and B, to zero: 
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Taking the partial derivative of eq. (3) with respect to a and B and 


setting it to zero yields: 





n _ n ey 
oA 2 _ pings Finx.+B MNBrYx*-p LY x%inx. = (4) 
ii ii tel 
OA _ no o c_®Q 
op mi B (a+1) ns (5) 


n 
n a 
na y xe 
B Ba = 1 
n esse (6) 
Ce ae 
i 
n 5 
B= (42 x] (7) 
i= | 
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Now replace the n in the 2nd term of eq.(4), with eq. (6): 


104 


n ee : “Ob as ay 
ze (5a2x7 jnb+ LInx +f InBxx*-B 2x*Inx, = 6 
=I i= i=l i=l 


and cancel with the 4th term: 


“ n es, 
ot Linx -B Yx*Inx =6 
= ] ] 1 


isl 


I 


From eq.(6), substitute 


5] 


Rearranging to set the right side equal to a: 


n 
Lx INx, . 
1 izt l 
i q "aa 
2X: 
i=l (0) 
= 1 
= n 
yx? Inx : 
i= | 1 
= = LInx 
be a i=l 


O72 


(8) 


APPENDIX F 


LIKELIHOOD DERIVATION FOR EXACT DATA WITH 


iF 
TYPE I AND TYPE I-V CENSORING 


For Exact data which has both Type I and Type I-v (Variable: Cause of 


failure not quantifiable in this experiment) censoring: 


Likelihood. =[[]f(x,.a.B)][][! -F(4,.0.B)][1-F(x,.0.8)] 


Exact Type I-v Type | 


: 1- F(x,,a,B) =ex 1 =) 


n = The total numberof fibers 
m =The total numberof fibers that failed 


m = The numberof fibers pamilioe © To) Ie 
m-m, = The number of fibers subjectto TypelI-v 


censoring 


a8 


Therefore, 


ce Ope 8 7 x exnf-Bxt} 
i=] 


eS) HlerHaN ‘ 


Let A=InL 


m m 
l =m 

=m Ina-m otnB+(a- 1) LInx - 8 Ee 
l=) i=] 


m X 
Bes (n= m)) (2) 
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To find the maximum likelihood of L, 


Taking the partial derivative of eq.(2) with respect to @ yields: 


a ot, al one Eee. ell 
= i inft 2D Iinx =8 2 x7inx.-fB iInBdx”| - 
= ist et ' ) ist | 


@e) 


pee "ing -B“InB y gi |- - m)() in| =0 (3) 
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Taking the partial derivative of eq.(2) with respect to B yields: 


m 
= ae eee, = ee, ox =) 
5B 77 m8 Sac eee se < +0” (n~m)x*=@ 
—| i 


1=mM 
] 


a+ 





Multiplying through by a 


m 


mt SX xX t(n —m)x=8 
i=l em, J 
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Now divide by my and rearrange to isolate B: 





To solve for a, first multiply eq. (3) by B% and rearrange: 


of M1, ol eal Te 
(3 se 1 UNE Linx |+InBd x*+InB X +(n -m)x@InB = 
(=a ) i= 1 1I=m > 
! 


m 
OL aA A OL 
Sx Inx + Dex Inx +(n- m)x‘Inx, 
ie 


: l 
i= i= 


Now gather the InB terms. 


mM 
| | 


4 a - m ee 
Bel eee ee ne Bee oe ie, ee In Y x's Ye eee 
i=] Re i=m , j 


mM 
| 


m 
AM A 
ne so 8 x Inx. +(n - m)x-INx 


l i : 
i=t =m 
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Divide both sides by m, and rearrange to bring & to the left hand side: 








m 
] m 
Aa A 
Batts, + » Ang, s¢0—mat in, 
i= I ; 





diel 
= a ml a a Aa oO is 
m, yx. + Lx. +(n —m)x? Inp 
ic1 ° ism, © J 
a 
In 1 s' (5) 
bee 2Inx. 
1 j=1 
Obtain from eq. (4), 
O a mi Aa 
mie = >x + > x. + (nm ~m)x” (6) 
] | | J 
| 1=M | 


and substitute for (a) in eq. (5): 


m B* [m B° [inp , 


a 


i 
InB + ay SU, 


m 
| m 
a A 
| poe a Ua oy x Inx +(n ~m)xINX , J 
] i= =m 


a! 


Cancel the InB terms and invert both sides: 


m 


] m 
Aa “A 
Dx IN. Se ne 
ict > 4) jem * 3 sow mn 
oy: = | ee 
= ~ | 
mp at > dee 


Finally, eliminate B by substituting eq. (6) for the denominator: 
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I 
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Ar An 
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Figure 3. Expected Rank of 008 strength set {D.} 
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Figure 4. Expected Rank of 019 strength set {D.} 
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Figure 5. Expected Rank of 008 strength sets {D,} and {D5} 
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Figure 6. Expected Rank of 019 strength sets {D.} and {D;.} 
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Figure 7a. Expected Rank data of merged 008 strength set {D.+D;.} 
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Figure 7b. Expected Rank data of merged 008 strength set {D.+D;5} 
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Figure 7h. Expected Rank data of merged 008 strength set {D.+D 5} 
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Figure 8e. Expected Rank data of merged 019 strength set {D.+D;<;} 
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Figure 8f. Expected Rank data of merged 019 strength set {D.+D;.} 
(cont.) 
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Figure 8g. Expected Rank data of merged 019 strength set {D.+D;.} 
(cont.) 
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Figure 8h. Expected Rank data of merged 019 strength set {D,+D 5} 
(cont.) 
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Figure 81. Expected Rank data of merged 019 strength set {D.+D;.} 
(cont.) 
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Figure 8). Expected Rank data of merged 019 strength set {D+D;.} 
(cont.) 
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Figure 8k. Expected Rank data of merged 019 strength set {D.+D};.} 
(cont.) 
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Figure 10. Likelihood surface and contours of 019 strength set {Dg} 
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5 % and 1 % confidence intervals for 008 strength set {Dg} 
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Figure 12. 5 % and 1 % confidence intervals for 019 strength set {D.} 
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Figure 13. 
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Figure 14. 5 % and 1 % confidence intervals for 019 strength sets {Dg} 
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Figure 15. Expected Rank of 008 strength sets {D, +D;.} and {D.} 
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Figure 16. Expected Rank of 019 strength sets {D, +D;,} and {D.} 
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Figure 17. 5 % and 1 % confidence intervals for 008 strength sets {Dg} 
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Figure 18. 5 % and 1 % confidence intervals for 019 strength sets {Dg} 
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Figure 19. Expected Rank of 008 strength sets {D, +D;,} and {Dij5] 
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Figure 20. Expected Rank of 008 strength sets {D, +D;.} and {D;9<} 
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Figure 21. Expected Rank of 008 strength sets {D, +D;,} and {D;z.} 
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Figure 22. Expected Rank of 019 strength sets {D, +D;,} and {D4j,} 
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Figure 23. Expected Rank of 019 strength sets {D, +D;.} and {D;9<} 
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Figure 24. Expected Rank of 019 strength sets {D, +D;.} and {D;z.} 
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Figure 25. Expected Rank of 008 strength sets {D, +D;.} and {D} 
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Figure 26. Expected Rank of 019 strength sets {D, +D;,} and {D} 
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Figure 27. Likelihood contours of 008 strength set {D} 
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Figure 28. Likelihood contours of 019 strength set {D} 
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Figure 29, Likelihood contours of 008 life set {D4} 
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Figure 30. Likelihood contours of 008 life set {D;>} 
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Figure 31. Likelihood contours of 008 life set {D;¢z} 
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Figure 32. Strength-Life Model for AS4-008 
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